MATRIX DIVISIBILITY SEQUENCES 

GUNTHER CORNELISSEN AND JONATHAN REYNOLDS 

Abstract. We show that many existing divisibility sequences can be seen as sequences of deter- 
minants of matrix divisibility sequences, which arise naturally as Jacobian matrices associated to 
C^ • groups of maps on affine spaces. 
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1. Introduction 

The most famous divisibility sequence is probably the Fibonacci sequence {Fn}n>i'- if fri di- 
^ . vides n, then Fm divides F„. This property is shared by other linear recurrent sequences |3|, such 

as any other Lucas sequence, and by higher degree recurrent sequences known as elliptic divisi- 
bility sequences Ill3l l23'l. Recent years have witnessed a revived and increasing interest in such 
sequences |[lTl[T2l[l4l[l5, 161 . alongside applications in cryptography 1 18, 211 and undecidability 
^ ' lUlllOl. In the current paper, we argue that - in a non-tautological way - behind each of these 

divisibility sequences lies hidden a naturally defined divisibility sequence of matrices, such that 
the given divisibility sequences occurs as the determinant of the sequence of matrices. 
CN ■ The plan for the paper is as follows: we shall first introduce the general notion of a matrix 

divisibility sequence indexed by a semigroup. Then we will see how a faithful representation of 
Q \ the semigroup by endomorphisms of an affine space gives rise to a matrix divisibility sequence, by 

CN ■ considering the Jacobian matrices of the endomorphisms. We will show how most of the commonly 

^ \ known divisibility sequences (mentioned briefly above) arise as determinants of matrix divisibility 

t^ ' sequences through interesting semigroups of endomorphisms of affine spaces, often associated to 

a representation of addition in an algebraic group. For example, Lucas sequences are associated to 
the 2x2 Borel group. We also construct the elliptic matrix divisibility sequence that underlies the 
usual elliptic divisibility sequences, and prove that it has primitive right matrix divisor classes. 
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2. Matrix divisibility sequences 

In this section, we introduce general matrix divisibility sequences over a ring S, indexed by a 
semigroup F, and we define primitive divisor classes of matrix divisibility sequences. 
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2.1. Definitions. Let S denote a commutative unital ring and Matrf(5) the ring of d x d matrices 
over 5. A (right) divisor class of a matrix M G Matd(S') is a coset GLcl{S) • M of M in the 
left quotient of Mat(i(S') by the invertible matrices GL^(5) over S. A matrix M is said to (right) 
divide a matrix N if there exists a matrix Q such that N = QM. If M (right) divides N, then any 
element of the divisor class of M also right divides A^. 

2.2. Example. An interesting special case of matrix divisibility is that of integer matrices (i.e., 
S = Z). In this case, the right divisor classes of a matrix M are in bijection with subgroups of the 
cokemel Z'^/M^Z" of left multiplication by the transpose M^ of M, cf. |[4l. 

We will only consider right division from now on, and hence frequently leave out "right" from 
the terminology. 

2.3. Definitions. Let (F, •) denote a (not necessarily commutative) semigroup. A divisibility se- 
quence of matrices over a commutative ring S, indexed by F, is a collection of matrices 

in Matrf(S'), such that if a right divides /3 in F, then Mq, right divides M^ in Matrf(S'). A primitive 
divisor class of a term Ma of such a sequence is a right divisor class of Ma that is not a right 
divisor class of any M^ for /3 a right divisor of a. 

If {Ma}aer is a matrix divisibility sequence, then 

{det(M„)}„6r 

is a divisibility sequence consisting of elements from the ring S. This is obvious from the multi- 
plicativity of the determinant. 

In general, divisibility of matrices is strictly stronger than divisibility of their determinants. For 
example, the matrices diag(l, 2) and diag(2, 1) are not right or left divisors of each other over the 
integers, although of course, their determinants are. Thus, it appears that the theory presented here 
is a strict superset of the existing one. Over a PID, divisibility of matrices is in general also stronger 
than divisibility of their individual elementary divisors in the Smith Normal Form. 

3. Matrix divisibility sequences arising from endomorphisms 

We produce a natural source of matrix divisibility sequences, as Jacobian matrices of endomor- 
phisms of affine space. 

3.1. Definitions. As before, let (F, •) denote a semigroup, S a commutative unital ring. Now let 

[•] : F -^ End(A^) : a ^ [a] 

denote a faithful representation of F into the group (under composition) of endomorphisms of 
affine d-space A^ over S (i.e., morphisms A| — )• A^). Let x G A'^(S") be a point in some 
ring extension S — )• S'. The matrix divisibility sequence associated to (F, [•]) is the sequence of 
Jacobians {Ja}a&T{x), with Ja and dx d matrix whose (i, j)-entry is given by 

{Ja)i,j ■.= d{[a]{x))i)/dxj. 
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The associated determinantal divisibility sequence is given by 

{det(jQ,)(a;)}Q,gr- 

3.2. Example. A trivial example: set F = (Z>o, •), and [n] : A^ —5- A.^ : x i— ;■ x". Then indeed 
[mn] = [m] o [n], and the associated (matrix) divisibility sequence is nx^^^. At x = 1, this is just 
the divisibility sequence of integers 1,2,3,.... 

The following facts are obvious, but they represent the basic idea in our definition: derivatives 
turn composition into multiplication. 

3.3. Proposition. A matrix divisibility sequence associated to (T, [•]) as before is indeed a matrix 
divisibility sequence: if a right divides (3 in T, then for any x G A"'(S"), the matrix Ja{x) right 
divides J/3{x) in the semigroup ofd x d-matrices Mat„(S'), and det{ Ja{x)) divides det(J/3(x)) in 
S. 

Proof. Write /3 = 7 • a in F. Then [/3] = [7] o [a] . The chain rule for the Jacobian matrix implies 
that for any x e A'^{S'), we have 

Jl3{x) = J^{[a]x) ■ Ja{x) 

in Matrf(S'). One can then simply take determinants of this identity. D 

3.4. Remark. We have included the case of a general semigroup F, instead of focussing on the 
(positive) integers as index set for the sequence, because some natural examples arise from elliptic 
curves with complex multiplication |22|, and even noncommutative semigroups occur naturally 
from supersingular elliptic curves over infinite fields of positive characteristic. 

3.5. Remark. A more general case would arise when one replaces affine space A*^ by an algebraic 
variety X. If [•] : F ^-)- End(X) is a representation, then one may consider the pullback of [a] to 
the tangent bundle 

d[a] : TX -^ TX, 
which then satisfies the chain rule 

d[ap]{x) = d[a\{/3x)od[(3]{x). 
Instead of taking a determinant, one may construct the highest exterior power 

det d[a] : /\^TX ^ /\'^ TX 

as automorphisms of the canonical bundle /\ TX. In general, however, there is no canonical 
choice for compatible coordinates in tangent spaces at different points (as there is on affine space), 
so that this does not lead to a "numerical" divisibility sequence. Therefore, we will not consider 
this more general setting here. 
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4. A construction of endomorphisms from algebraic groups 

A natural context for endomorphism representations is one that arises from the endomorphisms 
of a linear algebraic group, as follows. Let (G, +) denote an affine algebraic group over a field k, 
and let T C Endfc(G) denote a finitely generated semisubgroup of the algebraic group endomor- 
phisms of G. Fix an affine embedding of G into A"^. Choose generators 71, . . . , 7„ for the group, 
and fix an algebraic formula (74) for the action of the generators on the affine embedding, and fix 
an algebraic formula for the product and inverse in the group in the given embedding. Now define 
a representation [•] : T -> End(A'^) by [J2 aai]{xi, ...,Xd):=J2 O'iili) {xi,---, Xd), where Sa^ 
is computed using the given formulas for + and — in the group. 



4.1. Example. Example 13.21 fits into this framework, if we consider x 1— ;■ x™ as iterates of the 
multiplication map on the multiplicative group Gm- A more interesting example is the following: 

4.2. Example (Borel group and Lucas sequences). Consider the Borel group B of 2 x 2 matrices 
with the affine embedding 

and the multiplication formula 

{Xi,Yi,Zi)QiX2,Y2,Z2) := {XiX2,XiY2 + Y1Z2, Z1Z2), 

corresponding to the product of matrices, and a similar one for the inverse. Now, forn G N = F, 
consider the endomorphisms given by 

vn ^n 

[n]{X,Y,Z) = {X,Y,Z)Q---QiX,Y,Z) = iX^,Y ,Z"). 

^^ V ' A — Z 

n times 

The associated matrix divisibility sequences of Jacobians of [n] is 

/ nX"-i 

UX,Y,Z)=i YP{X,Z) ^^ YP{Z,X) 
\ nZ"-i 

with 

_ nX^-\X - Z) - {X^ - Z'-) 

^(^' ^) ^^zp ' 

and the associated determinant sequence is 

vn yn 



det(J„)(X,y,Z) = n^X'^'^Z 



X-Z ' 
an inocuous modification of the Lucas sequence for X and Z (and independent of Y). 

4.3. Example. Similarly, taking powers of matrices M G GL(2) leads to a determinantal divisi- 
bility sequence of the form 
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with a = tr(M) and /3 = tr^(M) — 4det(M). Here, when /3 = (i.e., the matrix has two identical 
eigenvalues), the formula should be understood in the limit as /3 — )■ 0, which gives n^(a/2)^('^~^). 
It could be interesting to consider the determinantal divisibility sequence of more exotic linear 
algebraic groups. 

One might wonder whether for Lucas sequences, one can do with one dimension less, but this 
is not even true for Mersenne sequences and general sets of endomorphisms of the affine line, as a 
simple integration proves: 

4.4. Proposition. A generalized Mersenne sequence {x'^ — l}n>i cannot occur as a matrix divis- 
ibility sequence associated to a set of endomorphisms of A}, i.e., in dimension d = 1. 



Proof. If so, then there are polynomials /„ such that 
By integration, we find that 



-I "Jn / N 



Jn\^j — ; Z X -\- Cfi 

n + 1 
for some constants c„, but then n^^ fn cannot be a representation, because it already fails to satisfy 
fmn = fm° fn (for example, deg{fm{fn{x))) = (m + l)(n + 1) / mn + 1 = deg(/„„)). D 

In connection with applications of divisibility sequences in logic, we record the following. Re- 
call that a subset X C Z'^ is called Diophantine if there exists an algebraic variety V defined over 
Z and a morphism -k: V ^ A'^ defined over Z, such that the image of the set of integral points of 
X is the given set: 7r(F(Z)) = X. 

4.5. Proposition. Suppose {M„}„gN is a matrix divisibility sequence that arises as above from 
an affine algebraic group G/Z, evaluated at a point P € G(Z). Then {Af„}„gN is a Diophantine 
subset ofTJ^ , and the associated determinant sequence {det(M„)} is a Diophantine subset ofL. 

Proof. By the David-Putnam-Robinson-Matijasevich theorem (see, e.g., |fl7l), Diophantine sets 
over Z are the same as recursively enumerable sets over Z. We prove that the set {M„} is recur- 
sively enumerable. The formula that expresses [n]x in algebraic terms, for a general point x G G, 
is computable in finite time on a Turing machine. The same holds for its Jacobian matrix. Hence 
also the values of the Jacobian matrices at P are computable in finite time. Now the set {Af„} 
can be enumerated by running through n. The same holds for the determinant sequence, since 
determinants are computable in finite time. D 

4.6. Remark. For a set of endomorphisms of a projective algebraic group (e.g., an abelian variety), 
one can use the general construction from Remark [331 One may also try to adapt the previous 
method from affine groups, by fixing an equation for addition in homogeneous coordinates and 
consider it on the affine cone over the group. (A particularly simple example of such a formula 
arises from the complete group law on the representation of an elliptic curve in Edwards form, cf. 
Il9ll . ||2< ) However, in general one will then only have a projective composition formula 

[aP]{P) = K^p{P)[a]mP), 
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for some functions Aq,^^ on G — from which the associated Jacobian matrix divisibility sequence 
will not in general be multiplicative, but rather satisfy 

If P G G{S') is a point whose F-orbit stays within a fixed affine chart, then it is possible to 
extend the previous method. 

Another approach to general divisibility sequences, based on generalized GCD's, is due to Sil- 
verman ll20l . For a further approach to (non-divisibility!) sequences in higher genus, see Cantor 
Il6l (where the r-th division polynomial is zero at a point P if and only if rP is in the theta-divisor 
— compare with [5 1 for another interpretation of these sequences). 

In the next section, we will use a slightly different method for elliptic curves, based on the theory 
of division polynomials. 

5. Matrix elliptic divisibility sequences: formal construction 

We will now show how elliptic divisibility sequences fit into the matrix divisibility picture, using 
division polynomials. Let E denote a cubic curve with projective equation 

Y^Z = X^ + AXZ'^ + BZ^. 

over the ring S = Z[A, B]. The non-singular points of E over any field containing S form a 
group. Multiplication on the non-singular points of this cubic curve can be expressed using classical 
division polynomials 

(j)nix) UJn{x,yy 



We refer to lIH, |[I1 and ll20l (ex. III.3.7) for the definition of these polynomials. Here, 0„ and 

only depend on x. We now consider the following map of affine 2-space 

[n]: A^^A^: (X, Z) ^ (z"Vn(f ), ^"'V^n(f ! 

The multiplicative property {mn)P = m{nP) translates to [mn] = [m] o [n] (compare |J7l, formula 
(5)), so that [•] indeed defines a faithful representation of F = N as a group of endomorphisms of 
affine 2-space A^. Hence the associated sequence of Jacobian matrices is a matrix divisibility 
sequence, and its determinant is a divisibility sequence in the usual sense. We now establish a 
formula for these sequences in terms of known division polynomials. For this, we first compute 
some partial derivatives: 

and 



dX{[n]-{X,Z)) 



,2 



XZ" -'<p'^{X/Z) + n'Z'' -'<Pn{X/Z) 



dZ 

= Z^'-\n'Zct>n{X/Z)-Xct>[XX/Z)). 
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Also 



dZ{[n]-{X,Z)) 

dX 
and 

dZ{[n]-{X,Z)) _ ^„2_2,^2 

dZ 
We conclude: 



z" -'^[Xx/z) 



Z" -\n^ZiJn{X/Z) - X<(X/Z)). 



5.1. Proposition. The sequence 

1 (X 7\- 7«^-2 ( Z(t>'^{X/Z) r?Z4>r,{XlZ) - X(t>'^{X/Z) 

"^ '^ - ^ V Zm'{X/Z) n^^l{X/Z)-X{^l)'{X/Z) 

is a matrix divisibility sequence, which we call a matrix elliptic divisibility sequence, with associ- 
ated so-called determinant elliptic divisibility sequence 

det{Jn){X,Z) = n^Z^^^"~^^W{(t)nM{X/Z), 

where W{f^ g) = f'g — fg' is the Wronskian determinant of two functions f, g, and ipn := tpn- ^ 

5.2. Remark. By Cassels' Theorem I in Q, the polynomial derivatives (jy^ix) and •ipn{x) have all 
their coefficients divisible by n; we conclude that the matrix Jn{X, Z) is divisible by the diagonal 
matrix diag(n,n). 

We can further simplify the Wronskian determinant in Proposition 15.11 as follows: by taking 
derivatives on both sides of 

/ / ^^ 4>n{x) 

x(n ■ [x, y)) = — 

i>n{x) 

we find that _ 

dx{n- {x,y)) _ H^((/)„, V'n.) 



dx ^2 

To use p-functions, we switch to classical Weierstrass form, by writing x = xi/36 and y = 
2/1 /432, so that (xi, yi) satisfies the Weierstrass equation in traditional form yf = 4xf — g2X — g^ 
for (72 = — 5184yl and g^ = —1866245, and we can write xi = p{z),yi = p'{z) for p the 
Weierstrass p-function of the corresponding lattice. Then 

dx(n ■ (x,y)) 1 dpinz) n , dz p'(nz) 

— — — p [nz)—- = n- 



dx 36 dx 36 dx p'{z 

which we further simplify to 

y{[n]{x,y)) 2n 1 ( mhn 

y i^2 V'n V "02 

so that we finally find 

5.3. Proposition. The determinant elliptic divisibility sequence from Proposition 15. i I equals 

detUX,Z) = n3z2("'-i)^(X/Z) = 2n3z2("'-i)^;!^(X/Z). D 

-02 •02 
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This result shows that every elliptic divisibility sequence occurs (up to passing to a field exten- 
sion to divide a given point by 2) as a determinant divisibility sequence. 

5.4. Remark. We have already seen how Lucas sequences arise from the 2x2 Borel group. Since 
all Lucas sequences also occur as elliptic divisibility sequence for singular cubics, we immediately 
find from the previous section that they, too, fit into this framework (|23|, Thm. 22.1). 

6. Matrix elliptic divisibility sequences: integral values and primitive divisors 

We now turn to the issue of actually substituting a rational point on the curve into these new 
sequences. 

6.1. Proposition. Suppose that P = {x, y) is a rational point of infinite order on an elliptic curve 
E/Q with chosen short Weierstrass equation with integral coefficients, and write x = a/b"^ in 
coprime integers a, b. The determinantal divisibility sequence 

det Jn{o-,b ) 

is integer valued, and has primitive prime divisors for n sufficiently large. 

Proof. First of all, we quote a result of Ayad (OJ) to the effect that if we write 



n 



with An, Bn coprime integers, then 



62" ^lia/b") = BlQn, 

where Qn is only divisible by primes p for which P is singular modulo p on the given model (so in 
particular, Q„ has only prime factors from the divisors of the discriminant A^; of the given curve). 
This means that 

det J„(a, 62) = n364("'-i)%i(a/62) = n'^ • Q'^, 

1p2 B2 

where Q^ has only prime divisors from A^;. 

Now Silverman has proven the elliptic analogue of Zsigmondy's theorem fl9l, implying that 
i?2n/^2 has a primitive prime divisor, say, p, for sufficiently large n (since P has infinite order in 
£^(Q)). We claim that p is coprime to n for n sufficiently large. Indeed, suppose p \ n. Since p is 
prime and primitive, P mod p has order 2n in E{Fp), so that by the Hasse-Weil bound 

2n <p+l + 2^ < n + 1 + 2^/n, 

leading to n < 6. Hence for n sufficiently large (n > 6, n large enough for Silverman's result 
to hold and for 2nP not to be an 5-integer, where S contains the primes dividing A^;), p is also 
primitive for det J„ (a, 6). D 

We finish this section by proving a matrix version of the existence of primitive divisors, based 
on the following general lemma: 

6.2. Lemma. Let {M„}„gN denote a matrix divisibility sequence in integral matrices Mn G 
Mat„(Z). If the associated determinantal divisibility sequence {det(M„)}„gN has primitive 
prime divisors, then the matrix divisibility sequence has primitive right divisor classes. 
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Proof. A nice way to organize the proof is by using the correspondence from Example 12.21 which 
implies that M„ has a primitive right divisor if and only if 7/^ /M^ T/- has a subgroup that is not 
in the image of any of the natural reduction maps !/■ jM^TJ^ — ;• TJ^ jM^l/- for any m \ n with 
m ^ n. But since we assume that det(M„) has a prime divisor p that doesn't divide any det(Mm) 
for any m \ n with m ^ n, p divides one of the elementary divisors of M„, but none of those of 
such Mm- This implies that the subgroup TJ^ jpTJ^ corresponding to p has non-trivial reduction, so 
corresponds to a primitive right divisor class. D 

6.3. Corollary (Elliptic matrix Zsigmondy theorem). Suppose that P = (x, y) is a rational point 
on an elliptic curve E/Q with chosen short Weierstrass equation with integral coefficients, and 
write X = a/b"^ in coprime integers a, b. There exists an integer N such that all the terms of 
a matrix elliptic divisibility sequence {Jn(a, 6^)}„gN with n > N have primitive right matrix 
divisors. 

Proof. This follows from the previous lemma since the associated determinantal sequence (cf. 
Proposition 15.11 ) has primitive prime divisors by Proposition 16.11 D 

6.4. Remark. One may also ask for 'converse theorems' in the following style: if the height of 
the entries of the matrices {M^} has a specific growth behaviour in n, does it follow (at least 
generically) that its determinant sequence has a 'related' growth behaviour? 

6.5. Remark. Linear and elliptic divisibility sequences satisfy recurrence relations (provided their 
terms are chosen with the right sign), so we ask: Is there a choice of representatives for the divisor 
classes corresponding to a Lucas or elliptic matrix divisibility sequence as in Proposition 15.11 such 
that these representative matrices themselves satisfy a polynomial recurrence relation (i.e., with co- 
efficients that do not depend on the index of the term of the sequence)? We have checked by direct 
computation that it is not the case that the "Borel" matrix sequence Jn{X, Y, Z) from Example 14.21 
satisfies a second order linear recurrence in matrices of the form 

Jn = A- Jn-l + B ■ Jn-2 

for matrices A = A{X, Y, Z) and B = B{X, Y, Z) independent ofn. One might argue that in the 
non-commutative ring of matrices, a second order linear recurrence should be of the form 

Jn = A- Jn~l ■ B + C ■ Jii-2 • D 

for matrices A, B, C, D independent of n, but we did not investigate this possibility any further. 
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